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ABSTRACT 

The spectral method of sea and swell wave forecasting yields more 
meaningful results than do other forecasting methods. To shorten the time 
required to make such a forecast the authors used high-speed computer 
methods and the Pierson-Neumann-James spectra to make spectral forecasts 
and display the results in graphic form for several synoptic weather models 
and a wide range of wind and fetch speeds. Forecasts were completed for 
a stationary fetch model, a fetch moving to leeward and a fetch moving to 
windward.— 

Computations were made with the assumptions that a rectangular genera- 
ting area could be delineated; that a spatially-uniform mean wind existed; 
that this wind did not increase or decrease until the time of forecast; 
and finally that if the fetch were moving, its movement was steady. 

The following parameter limits were established: 

a. wind speeds, 15 to 55 knots; 


b. speed of the generating area, 15 knots to and including 
the wind speed; 


c. forecast times, 6 to 72 hours. 
Results of the forecasts are plotted on approximately 250 graphs which 
show the geographic distribution of spectral sea or swell components 


relative to the fetch at forecast times. 
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1. Introduction 


It is universally understood that a need exists for a quick, yet 


accurate method for forecasting the wave conditions on the surface of 


the oceans of our planet. However it must be strongly emphasized that 


such a method should be spectral. It is important that the range of 


periods as well as heights of the waves be forecast. Some examples of 


where knowledge of wave periods is important are as follows: 


(1) 
(2) 
(3) 
& 


the emergency ditching of aircraft at sea, 
seaplane and aircraft carrier flight operations, 
refraction and shoaling of waves as they approach a shoreline, 


submarine surfacing and submergence. 


This need for rapid spectral wave forecasting prompted the authors of 


this paper to perform the study and work herein discussed. 
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2. Background 


Generally two catagories of forecast methods are available, a spectral 
method and a non-spectral method. The spectral forecast has the advantage 
of determining the range of periods with important energy but is a very 
tedious, time-consuming process. Furthermore, the technique requires that 
the forecaster be trained quite extensively. The non-spectral method has 
an advantage in that a forecast can be made rather rapidly and by a person 
with minimal training. However, it does not tell the user of the forecast 
the range of periods in which the wave energy lies. 

Therefore, in order to permit accurate and more meaningful spectral 
forecasts by relatively untrained persons, high-speed computer methods 
were used to make spectral forecasts for as many synoptic weather models 
as time would permit and to display the results of the computer output in 


graphic form. 





3. Basis of the Forecast Method 

The Pierson-Neumann-James forecast method as outlined in H. O. Pub. 
No. 603 [3] was chosen as the basis for the wave forecasts. In order to 
avoid repeating a great volume of theory and methods, it will be assumed 
that the reader has a working knowledge of this publication. It is recog- 
nized that Pierson and Moskowitz [1 and 2] are developing new spectra 
which will improve upon the PNJ spectra now in widespread use. Therefore, 
the computer programs were written in such a manner as to allow utiliza- 
tion of these new spectra or others as they are published. 

In making a forecast of sea and swell conditions, one must know at 
what location within the generating area and how long after the wind began 
to blow that waves of a given period are first formed. This information 
is gleaned from the Co-Cumulative Spectra. Then one must know the velocity 
with which the waves are propagated. Each spectral component moves at a 


its respective group velocity, the direction being accounted for by an \ 


jf 


f 


angular spreading function. When these characteristics have been determined 


ee 


re 


the distribution of wave periods over any geographical area can be shown 
graphically. 
For each wind speed and fetch description considered, this distribu- 


tion is shown for various times since the wind started blowing. 


liane a ait es 
—e an ee a <> 
a) ie = f= « 


=i © aa © 
=> SS SS 


=_> ae 





4. Synoptic Models Completed 


To date, three models have been completed: 


(1) a stationary fetch model which may be associated with 
a stationary cyclone, tiough or ridge; 


(2) a model with a fetch moving to leeward, as associated with 
the southern quadrant of a cyclone moving eastward; 


(3) a model with a fetch moving to windward as associated 
with the northern quadrant of a cyclone moving eastward. 


In the following sections the variables and equations associated with 


each model are discussed. 


A discussion of the computer programming and procedures is found in 


Appendix A. 





5. Stationary Fetch Model 
5.1. FOPRIM® » 


With this case a description of wave conditions within the generating 
area (SEA) must be given. In order to generate waves of a given period, 


the wind must blow over a minimum length of fetch for a sufficient time. 


This minimum distance, measured from L e_fetch, is 


given the name\FOPRIMEJ The values for each period are determined from 


the CCS, 


FOPRIME = F '(V,t,) Ci) 


where V is the wind velocity and ty is the duration time. Graphs for 
specified velocities and various times of observation show plots of the 
highest period present versus FOPRIME. These graphs are also called 
FOPRIME. The use of these and all other graphs is explained in Section 10. 
5.2 SWELL 

A description of the waves emanating from the generating area (SWELL) 


in the general direction of the wind also must be given. As soon as a 


Spectral component with period T, is formed, it_moves—with—a-group vetocity, 


gee Dt ET 


¥. = 1.51 T, (V,t4) (2) 





where T., the _period_of the longest wave, is a function of wind velocity 

] ma staat tt entnennettmeieetiiminntemainend! 
and duration. Assuming that all swell waves are propagated from the lee- 
ward edge of the fetch, their maximum distance from this edge at some obser- 


vation time t is * 
ob 


ne ee 


SWELL = (top - tyV.: (3) 
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Graphs of spectral component versus SWELL for various observation times 


are also called SWELL. 





6. Moving (leeward) Fetch Model 
6.1 FOPRIME 

For this synoptic model the description of the SEA conditions within 
the windward and leeward portions of the fetch must be given. Since the 
fetch is moving, it is assumed that the windward portion of the fetch is 
completely filled in with waves that have group velocities less than the 
fetch velocity. FOPRIME for each spectral component with group velocity 
greater than the fetch velocity is obtained from the CCS, The graphs 


showing the distribution of spectral components versus FOPRIME are called 


FOPRIME. 
6.2 SUMDEL 


Due to the_movement of the fetch into an area where it—is—assumed— 


calm seas exist, the effective duration of the wind blowing over the lee- 
feat a Pe ; Saar 
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ward portion of the fetch is less than for a stationary fetch. Therefore, 
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the distance behind the leading edge of the fetch that a spectral com- 
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ponent of period, T., can first be formed during an increment of duration 
— cL ae eS ae ll 
time At. is 
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aD, = (Vp-1.515 ae (4) 
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where V; is the fetch velocity. Each spectral component requires a specific 
eneration time; a summation of the duration increments At, determines this 


N SS a 
‘generation time; t, =, 2,4t, . The summation of AD's determines the total 


lag of that component behind the leading edge of the fetch; this distance 
is called SUMDEL. 


N 
SUMDEL =,2, AD, . (5) 


If a spectral component generated at a particular time t 
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velocity greater than the fetch velocity, then at time i, this component 
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has advanced’ relative to the leading edge of the fetch an amount, | 
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where DIST is the travel distance of the component since it was generated 
as a alr pt IAA ET a I GE Atte A Ee 








Fad Ss SS 


ae Srenteneneecns ae a : ae aE 


) Ve is the distance the fetch advances during this same period 


DIST = (t ty) Ue (6) 


Hence SUMDEL is modified by this relative travel of the high-speed spectral 
component. The resulting distance that a specific spectral component lags 


behind the leading edge of the moving fetch at t | is called SWEDEL. 


b 


SWEDEL = SUMDEL - DIST + (t pty) Vee (7) 


6.3 SWELL 


The distribution of the swell to leeward of the fetch is similar to that 





in the stationary fetch model, but since the fetch is moving it becomes 





necessary to consider the distance traveled by the leading edge since the 
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generation of a specific spectral component, (cme ta) Ve» and SUMDEL, the 
E Prine: a | 


distance behind the leading edge that a specific component is formed . Thus , 


in this case, 


say US = (t .-t.) V, = Soin tae V_ ~ SUMDEL. (8) 
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Graphs showing the range of periods of the spectral components versus SWELL 
for various observation times are called SWELL. 
6.4 LAG 

As the fetch moves to leeward at a specific velocity, spectral compon- 
ents with group velocities lower than the fetch velocity fall out of the 
windward portion of the fetch, assuming as before that the windward part of 
the fetch is completely filled in by these spectral components. The distance 
that a wave component with period T, lags behind the trailing edge of the 


fetch at a is termed LAG, where 


b 
LAG = (t opty? Ve - (t pty) Vi» (9) 





or by rearranging terms 


LAG = (t ptg (V_ 7 V2 : (10) 


Graphs showing the distribution of spectral components versus LAG are 


called LAG. 





7. Moving (windward) Fetch Model 

For this situation only two of the areas previously mentioned need 
‘describing, the area within the fetch and the area behind the trailing 
edge where swell is found. 
7.1 SUMDEL 

The distance behind the leading edge of the fetch that a specific 
spectral component is first formed is called SUMDEL as defined by (5). 
Since the fetch movement and wind direction are opposite, (4) is modified 
to 


AD. = (V; + 1.515 T,) At. (11) 


The graphs describing the generating area in this case are called SUMDEL. 
7.2 LAG 

Immediately upon generation, waves of a given spectral component move 
out of the generating area. Their distance behind the trailing edge is 
called LAG. Because of the opposite directions of movement of fetch and 
wind, equation (10) becomes 


LAG = (t n7tg) Ve + V.)- (12) 


The graphs depicting the location of spectral components behind the fetch 


at a specific observation time are called SWELL. 
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S. Definitions Used in the Forecast Method 
To make clear the notations used in subsequent sections, the follow- 
ing terms are defined. Most terms are the same as those used in [3]. 
However, some of the terms are unique to this paper. 
Fetch - That limited or bounded area of the ocean where it is 
assumed a mean steady wind is blowing and generating waves. 


The area is also assumed to be rectangular with one axis 
along the wind path. Also referred to as a generating 


area, 

V Sea - Those waves found within the generating area. 

V swell - Those waves which have moved out of a generating area, 
Vv R - The distance from the center of the leeward edge of a 


generating area to the point P at which a forecast will 
V Ke be made. 


L - The distance from the center of the windward edge of a 

O ; : 
moving (leeward) fetch to P, measured in a direction 

V 7 opposite to fetch motion. 

Ty - The highest wave period in either a sea or swell. Ty 
is a function of wind speed and duration. 

Tey - The highest wave period in either a sea or swell where 
a relatively short fetch limits the growth of waves. 
Refer to section 13.1d. 

qT, - The lowest period in either a sea or swell. 

Ey - The energy in that part of the spectrum between the periods 
zero and Ty: 

EF. - The energy in that part of the spectrum between the periods 
zero and T, . 

AE - The energy associated with that portion of the wave spectrum 
with periods from Ty to Ty AE = Ey 2 EF. 


8; TA - The angles between the direction of the mean wind and the 
lines drawn from Pto the center of each side of the fetch. 


S - The angular spreading factor, determined by 0, and e, from 
the Angular Spreading Factor Tabie, explained in Section ll. 


AE - The total increment of energy arriving at P at a specific time. 
This energy is a function of wind speed, fetch speed, distance 
from the fetch and the relative direction of P from the fetch. 


AE, = AE(S,) = (Ey = E,) (Sp) 
10 





t b - The time of the pan si measured from the time the wind 
i began to blow (e. g. if the wind begins at 010000Z and a 
forecast at P is we be made for 020600Z then t i is 30 hours). 
Also referred to as observation time. 
V - Mean wind speed in knots. 


- Mean fetch speed in knots. 


ri 





9. Parameter Limits 

The range of wind velocities considered for the computer computations 
was from 15 to 55 knots. Below 15 knots, the wave energy is considered 
insignificant. The upper limit of 55 knots is controlled by the fact 
that the PNJ spectra are given only up to 56 knots. Since facsimile 
weather maps show wind speeds in increments of five knots, this increment 
value was used for the computations. 

The range of fetch velocities used was from 15 knots up to the existing 
wind velocity. It was considered unlikely that a delineated generating 
area would move faster than the wind velocity within it. For fetch velocities 
less than 15 knots, it was considered that the stationary fetch graphs could 
be used, assuming some mean location of the fetch. 

Graphs to forecast sea and swell conditions were drawn for observation 
times from 6 to 72 hours. The upper limit was arbitrarily set by the authors. 
The observation time was incremented by six hours for convenience. Any other 
times can be found by interpolation. 

For the moving-fetch models, a fetch-length limitation of 1600 
nautical miles was established for reasonable graph scaling. This limita- 
tion seems reasonable since it is doubtful that a uniform wind field larger 
than this would be found. This fetch limitation is especially noticeable 
in limiting the formation of low-frequency waves in the Moving (windward) 


Petch.Model. 





10. Use of Graphs 


Each of the graphs discussed has 12 separate or superimposed curves 
showing the range of wave periods at geographic locations for observation 
times from tub 6 hours to tb 72 hours at six-hour intervals. Fig. 1 
shows a time series of graphs approximating a situation where a generating 
area with a mean wind velocity of 35 knots (V = 35) moves across the ocean 
with a fetch velocity of 20 knots (V >= 20). The final graph in each column 
is the composite picture where curves for all times are superimposed. The 
shaded areas in each case show the range of periods present at all locations. 
All distances measured on the graphs are relative to the windward or leeward 
edge of the fetch. 

To make forecasts for times other than those of the observation time 


curves, the forecaster may interpolate between curves. 
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Figure 1. Time Series of Graphs Showing Distribution 
of Spectral Wave Components Relative to 
a Moving (leeward) Fetch. 
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11. Use of the Angular Spreading Table 


The angles 8., and 9, defined in Section 8 determine the reduction of 
energy due to the angular spreading the waves undergo after leaving the 
generating area. The angular spreading factor nomograph in [3] was used 
to formulate Table 1, showing the angular spreading factor So as a function 
of the two angles in ten-degree increments. It is felt that measuring these 


angles to the nearest ten degrees gives sufficiently accurate results. 


15 
















Sere P| fel pe bela be peer 
[29] - [21 [02 [07 | | 20 20 | 40] 50 2070 20 |20 2 [7 90 [to 
a red eeteeeeee et elec 








= 
[fal 

Sa ot 
wg gS MS aed al ED Sa 
ire ed fal EES Eo te a Re EE EE EE 
peyee Hem [2a false 0 fs | 0] 40 20] 40 |n0) 28 | 2 
Ta Fed ind Fob el El 1 ES EEE 
peepee fer in| fe Jeet oh 
| 70 for | [92] fo | tf se = 0 | 10 [2 | 27 
120] 2070 [77 | 79 | s0} 50] 50] 40| 20] 20 [10 | - [os [13 [17 [10 [20 
[40] 2° |e? Jos [22 | 76 [50 | s0] 40] 20].20| 10 [00] - [os [oo |m |e 
-sa} 2 [2 [oo eee aise 5. outer 





mailer eet 
2020 | 70} 20} 20] 0} 20] 20] se Jor] os] or] — 
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12. Use of the Energy Table 
For a rapid forecast, measuring the range of wave periods to the 
nearest whole second at any given location seems sufficiently accurate. 
Hence the simplified Energy Table (Table 2) gives the energy associated 
with a spectrum of waves with periods from zero to T, listed on the left 
vertical edge of the “sble. Enter this table with V, and Ty oF T to 
determine Ey or EF respectively. For more accuracy, the period to the 


nearest tenth may be read from the graphs and the energy value from Table 


2 interpolated. These extracted values have units, feet . 
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13. Forecasting Procedures 


The following is a concise, step-by-step procedure to be followed in 


making a forecast. 


13.1 Preliminary steps 


a. 


Draw the boundaries of the generating area on a synoptic weather 
chart; determine the mean wind, V, within the area and if the 


fetch is moving determine a mean fetch speed, V 


f° 


Locate the point of interest, P, where the forecast is to be made. 


Determine into which of the synoptic models the forecast falls: 


i 
Z. 


3: 


Stationary Fetch Model (Section 13.2) 
Moving (leeward) Fetch Model (Section 13.3) 


Moving (windward) Fetch Model (Section 13.4) 


13.2 Stationary Fetch Model 


A sea or swell forecast is made for P depending on the location rela- 


tive to the fetch as shown in Fig. 2. 


a. 


ie 


Ia 





Figure 2. Sea and Swell Areas Associated with a 


Stacionary Fetch 


Sea Forecast 


Measure the distance from P to the windward edge of the fetch. 


Enter the FOPRIME graph for the appropriate windspeed with 


this distance and t to obtain T... 
ob U 


Assume Ty is zero; therefore, ES is zero. 


From Table 2 determine Ey a function of Ty and V. 
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OF 


Assume that within the generating area the angular spreading 
factor Sg is unity; therefore, OE equals Ey: 

For characteristic wave height computation, refer to the 
formulae in Table 3. 


b. Swell Forecast 


(If P lies outside the sector outlined in Fig. 2, swell will be 
considered insignificant.) [See Note] 


l. 


dee 


5. 


Determine as outlined below if the length of the fetch will 
limit the growth of waves (e.g., a 30-knot wind with a 
duration of 24 hours can produce waves with periods up to 
16.7 seconds if the fetch length is at least 280 nautical 
miles. If the fetch length is only 200 NM, the highest 
period of wave which can be generated is 11.7 seconds). 
Enter the FOPRIME graph for the given V, V, and t b with 
fetch length to determine the highest period of waves 
generated, If this period is less than the highest period 
on that t_, curve, then it is called T,,, and replaces T 

in the remainder of the swell forecast. (See Example Fore- 
cast No. 1.) 


Measure Ro: 


Enter SWELL graph for appropriate V and t _, and obtain T 
ob U 
(or Tay and TL. 


From Table 2 determine Ey and EL: 


Draw lines from P to the centers of the sides of the generat- 
ing area and measure angles 9, and 8, as shown in Fig. 6. 
Angles are negative if measured clockwise from the wind 
direction, and positive if measured counterclockwise. 


Enter Table 1 with 9, and 9), and determine Sy: 
Calculate AE, = (Ey - E,) (Sy) 


Use formulae in Table 3 to determine wave heights. 


13.3 Moving (leeward) Fetch Model 


Fig. 3 is an aid in determining the type of waves to be forecast. Prior 


to any other steps the fetch must be tracked ahead to its location at the 


time of forecast, t 


ob’ 


Note: When a long generating area exists, swell may emanate from the sides 


of the fetch. 


In such cases, the 45°-sectors shown on Fig. 2 may be moved 


toward the windward edge of the fetch an amount deemed necessary by the 


forecaster. 
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Figure 3. Sea and Swell Areas Associated with a 


Moving (leeward) Fetch. 
Sea Forecast 


(Since the graphs are drawn to cover a wide range of conditions, 


the highest period of sea waves may be determined from either the FOPRIME 


or SUMDEL graph. Generally, for slow fetch speeds T,, will be determined 


U 


from the FOPRIME graph and for high fetch speeds T,, will be determined 


U 


from the SUMDEL graph; therefore, both graphs should be investigated. ) 


1. Measure the distance from P to the trailing edge of the fetch. 


2. Enter the appropriate FOPRIME graph with this distance and 
Cob to obtain Ty: 


3. Measure the distance from P to the leading edge of the fetch. 


4, Enter the appropriate SUMDEL graph with this distance and 
t_, to obtain T,. If two separate bands of wave periods 
fe) 
are in evidence, assume all periods between zero and Ty 
are present. 


5. Using the smaller of the two values of T 


2 and 4, enter Table 2 and obtain Ey: 


6. Assume T, is zero; therefore, EF is zero. 


7. Assume S4 = 1; therefore OE, = —E 


8. Use the formulae in Table 3 to determine wave heights. 


y obtained in steps 


i 


Swell Forecast for Area Leeward of the Fetch 
(Swell is present ahead of the generating area only in the sector 


shown on Fig. 3 and only if the fetch is moving at speeds of 20 knots or 


less.) 
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1. As in the Stationary Fetch Model forecast (Section 13.2b), 
it must be determined if the fetch length limits the gener- 
ation of waves. With fetch length enter the FOPRIME and SUMDEL 
graphs for the given V, V, and t , to determine Tey from each 
graph and select the lower of these values. 


2. Measure Roe 


Ww 


Enter the appropriate SWELL graph with R_ andt b to obtain 
Tei Tee andaT.. ‘ ‘ 

U FU L 

» From Table 2 determine Ey and Eye 

oi ) 


4 

5. Determine 60 o,, and S.., 
6. Calculate AE» = (Ey - ED (So). 
i 


. Use formulae in Table 3 to determine wave heights. 


c. Swell Forecast for Area Windward of the Fetch 


U 

. Since it is difficult to determine an angular spreading 

factor for swell waves traveling slower than the fetch, 

assume S, = 1.0 and hence the swell heights are the max- 


imum possible for this area. Hence, AE, = Ey: 


1. Measure L.- 

2. Enter the appropriate LAG graph with L, and t ob to obtain Ty 
3. Assume T, is zero; therefore, Ey is zero. 

4. From Table 2 determine E,.. 

5 


6. Refer to Table 3 for height computations. 
13.4 Moving (windward) Fetch Model 
Fig. 4 describes the location of swell relative to the generating area. 
Prior to any forecasting steps the fetch must be tracked ahead to its loca- 


tion at Cob: 





Figure 4. Sea and Swell Areas Associated with a 
Moving (windward) Fetch. 


Z2 





a. Sea Forecast 


i 


6. 


Measure the distance from P to the leading edge of the 
fetch. 


Enter the appropriate SUMDEL graph with this distance 


and Cob to determine Ty: 

From Table 2 determine Ey: 

Assume To is zero; therefore EF is zero. 
Assume S 6 = 1.0; therefore, AE = Ey 


Refer to Table 3 for height computations. 


b. Swell Forecast 


hes 


6. 


= 


Determine if the fetch length limits the growth of waves 


by entering the appropriate SUMDEL graph with fetch 
length and t ob as in Section 13.2b. 
Measure R.. 

Oo 


Enter the appropriate SWELL graph with Ro and t Ob to 
obtain Ty (or Tay? and TL. 


From Table 2 determine Ey and ES - 


Measure 9. and 9), and obtain So from Table 1. 


Calculate AE, = (Ey - E) (S,). 


Refer to Table 3 for height computations. 


13.5 Multiple Storms 


Should it appear that the forecast location might be under the 


influence of more than one storm, the incrementsof energy, AE 9 » at P 


from each storm should be calculated individually and then summed. 


characteristics are computed from this summed energy increment. 
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14. Wave Heights 


The following formulae used most often are the same as those appear- 
ing in [3]; they have been modified only to agree with the symbols used 
in this paper. For determining range of heights, exceptionally high waves, 


and amplitude information, refer to [3], Chapter 1. 


Table 3. Wave Heights. 


Most frequent height . .....e « «e « « See = 1.41VAE, 
Average wave height. . . . « « « o e « « ae = 1.77/4E, 


Significant wave height (highest 1/3). . H 


Average height of highest 1/10..... Hy 10 = 3.60/AE, 
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15. Forecast Worksheet 

In order to systematize the described processes, a worksheet (Fig. 5) 
was devised. The blanks provided are filled with information extracted | 
from the graphs and tables of the system. It is recommended that any unit 
using the forecast method reproduce this worksheet. The completed work- 


sheet may be kept on file to provide a permanent record. 
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16. Example Forecasts 


Example forecasts, made for three synoptic fetch models, are organized 
in the format of Section 13 (Forecasting Procedures) with correspondingly 
numbered steps. Each step is amplified by reference to a figure or table 


from which the information is obtained. 


Ly 





16.1 Example Forecast No. 1: Stationary Fetch Model 








K— 460 wm—| 


Figure 6. Geographic Positions of P 


and P, Relative to a 
Stationary Fetch. 


1 


Fig. 6 portrays a fetch located in one geographic position for a 
period of at least 60 hours. From the synoptic and prognostic weather 
charts and Fig. 6 the following preliminary information is obtained for 
the forecast: 

1. the positions of Py and P53 
2. V =55 knots beginning at 1921002; 
3. a Stationary Fetch Model. 
Proceeding with the steps as outlined in Section 13.2: 


a. Sea Forecast for PS at 2200002; ae = 60 hrs. 


1. Distance from the windward edge to Pi = 275 NM 


2. Ty = 12 sec. (FOPRIME graph; Fig. 7) 
3 TT, = 0; EL = 0 

4 Ey = 64.0 (Table 2) 

Die AE, = 64.0 


So TA ye 2.83/64 = 2.83(8) = 22.6 ft. (Table 3) 


b. Swell Forecast for P, at 220000Z; t , = 60 hrs. 


2 ob 
1. No fetch limitation on highest period generated. (See Note 1) 





INO 
nw 
Il 


7OONM (Fig. 6) 


Lo 
J 
! 


= 19; Tt = 12 (SWELL graph: Fig. 8) 
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4. EB = 128; E, = 64 (Table 2) 


U 
Be. = +20°; eo, = +40° (Fig. 6) 
6. S = .18 (Table 1) 


Te AE = (E,,-E,) (Sy) = (128-64) (¢.18) = 64(.18) = 11.5 
8. H/37 2.63/11 .59— 2.83(3.4) = 9.6 ft. Gable 3) 
Note 1: If the fetch length were only 400 NM the building of waves 


would be limited so that T,,, would be 16 seconds (FOPRIME graph: Fig. 7). 


FU 


Then T,. = T = 7.6 ft. 


U FU = 16 seconds at P, and H 


1/3 
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16.2 Example Forecast No. 2: Movin leeward) Fetch Model 
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Figure. 9. Geographic Positions of P,_, Py and P, Relative 


to a Moving (leeward) Fetch. : 


Fig. 9 portrays the movement of a fetch from the time the wind began 


blowing to the time of the forecast, showing the three points for which 


forecasts are to be made following the procedure as outlined in Section 


13.3. The preliminary steps show: 


Me 


Forecast for P 


i. Py and P3> 


= 15 knots since 161800Z; 


the positions of P 


V = 40 knots; Ve 


a Moving (leeward) Fetch Model. 


at 181800Z: t = 48 hrs. 


l ob 


Distance from P, to trailing edge = 450 NM (Fig. 


1 


Ty = 14 sec. (FOPRIME graph: Fig. 10) 


Distance from P, to leading edge = 150 NM (Fig. 


i 
a not affected by SUMDEL in this situation. 


Ey = 130 (Table 2) 


AE, = E,, = 130 


1 


b. Swell Forecast for P, at 1818002: t ob = 48 hrs. 


ite 


2 


H hee 2.83/130 = 2.83¢11.4) = 32.2 ft. (Table 3) 


9) 


9) 


To determine if the length of the fetch limits the highest 
wave period which can be generated enter both the FOPRIME graph 
(Fig. 10) and the SUMDEL graph (Fig. 11) with fetch length 


and obtain: 
Bz 





Try = 18 sec. (FOPRIME graph) 
Tey = 22 sec. (SUMDEL graph) 
Using the smaller of the two, the highest possible wave period 
generated, Ty = 18 sec. 
U 
2, R, = 100 NM (Fig. 9) 
a Ty = Tey = 18 sec. (Step 1); qt = 14 sec. (SWELL graph: 
Fig. 12 ) 
4. Ey = 205; EY = 130. (Table 2) 
eG 430°; o, = 230° : S, = .6 (Fig. 9 and Table 1) 
6. AE = (205-130) (.6) = 45.0 
Les Ha /3 = 2.83/45 = 19 ft. (Table 3) 
c. Swell Forecast for P3 at 1818002: t Ob = 48 hrs. 
de, L, = 400 NM (Fig. 9) 
one Ty = 6 sec. (IAG graph: Fig. 13) 
oe Ty = 0; EL = 0 
4. Ey = 7.8 (Table 2) 
a, Eg = Ey = 7.8 
6. H/3 =2.83/7/7.8 = 7.9 £t? (Table 3) 
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16.3 Example Forecast No. 3: Moving (windward) Fetch Model. 
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Figtre 14. Geographic Positions of P, and P, Relative 
to a Moving (windward) Peeh 


° Pig. 14 portrays the movement of a generating area during the times 
indicated. From synoptic and prognostic charts the following preliminary 
information is obtained for the forecast; 

1. the positions of Py and Py 3 


2. V = 30 knots; V, = 20 knots since 0212002; 


f 
3. aMoving (windward) Fetch Model. 


Proceeding with steps outlined in Section 13.4: 


a. Sea Forecast for P) at 0400002; t ,. = 36 hrs. 
1. Distance from Py to the leading edge = 400 NM 
Ze T = 9 seconds (SUMDEL graph : Fig. 15) 
3 Ee 21.0 (Table 2) | 
4. T = 0; Ey =O, 
oe AE, = 21.0 


6. H/3 = 2,.83/21.0 = 2.83(4.6) = 13 ft. (Table 3) 


b. Swell Forecast for P, at 0400002; t ob = 36 hrs. 
1 Tey = 10 seconds (SUMDEL graph: Fig. 15) 
2 R, = 450 NM (Fig. 14) 
3. T, = TL, = 10 seconds (SUMDEL graph: Fig. 15 and SWELL graph: 
U FU ; 
Fig. 16) 
T = 0 


38 





eo 
ul 


29.0; (Table 2) 


. 
Eee 

0,20, 0, = +40°, S, = .38 (Fig. 14 and Table 1) 
AE, = (29) (.38) = 11.0 


Hi /3 = 2,83/11.0 


973 it. “(lable 3) 
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17. Further Investigations 


Following is a list of subjects which the authors hope to study further 
as time and computer facilities permit: 
1. completion of models where the wind is blowing at some angle 
with respect to the direction of fetch movement other than 
zero or 180 ; 


2. extending the duration beyond 72 hours; 


3. completion of a model for a dying sea and one with variable 
winds ; 


4. programming in SCRAP computer language for adaptation to U. S. 
Navy Fleet Numerical Weather Facility forecasting. 
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18. Graphs 
Included in this section are the graphs of LAG, FOPRIME, SUMDEL 


and SWELL, divided into sub-sections: 
18.1 Stationary Fetch Model; pages 44-62 
18.2 Moving (leeward) Fetch Model; pages 63-212 
18.3 Moving (windward) Fetch Model; pages 213-303 
The graphs within the sub-sections are arranged by increasing wind speeds 


and associated increasing fetch speeds. 
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18.1 Stationary Fetch Model Graphs 
INDEX OF GRAPHS 
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APPENDIX A 

iL. Programming, Computing and Graphing Procedures 

Fortran computer language utilizing the Control Data Corporation 
1604 computer was used exclusively. 

A brief outline of the procedure employed follows. 
1.1 The Co-Cumulative Spectrum arrays 

In order to have a table of values of maximum period and minimum required 
fetch as functions of wind velocity and generation time, the Co-Cumulative 
Spectra of Pierson, Neumann and James were read into the program as fifty- 
five by seventy-two arrays. That is, periods and Fo' were read into the 
computer memory for each five knots of wind from fifteen to fifty-five 
knots and every two hours duration time from two to seventy-two hours. 
Fig. A-l is a schematic representation of this array for period, T(V,t4)3 
where (V) is the wind velocity and ta is the duration. A similar array 


for Fo' was also read into the computer memory. 


Coisieg: A so eo ieee 
*30.,2 T50,4 NG Pee 2 ee MRM Te 
B55 "5 r,s 
55,2 ine Maer moe ee te 8 ee) 


Figure A-l. Wave Period Array. 
1.2 Computations for LAG, SWELL, FOPRIME, SUMDEL or SWEDEL 
Each model required the writing of a separate program utilizing the 
equations outlined in Sections 5, 6, and 7. Since the programs are unique 
to this thesis and would be of interest only to those readers wishing to 


pursue further study of these methods, they are not included. They may be 
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obtained by contacting the authors through the thesis advisor at the U. S. 
Naval Postgraduate School, Monterey, California. 
1.3 All condition print-out 

In order to check the validity of the computations, the programs for 
each of the three situations were run for the nine values of fetch velocity, 
fifteen to fifty-five knots, with the wind velocity ranging from the speci- 
fied fetch velocity up to and including fifty-five knots. When the validity 
of the computations was verified by random hand analysis, these computed 
values were then used to set the scaling factor for the proposed graphs of 
the four variables (LAG, SWELL, FOPRIME and SUMDEL). 
1.4 Modification of the Program 

Two modifications were necessary in order to use the all-purpose program 
in the computations required and its adaptation to a graphical output. The 
first step was the addition of the "call graph" statements. These state- 
ments called into the program the subroutine GRAPH. 
1.5 The Production Run 

For each situation the programs were run for the nine fetch velocities. 
Each run resulted in a tabulated output (an example of which is enclosed 
as Fig. A-2) and a graph-output tape. The graph-output tape was then run 


on the CDC 160 computer in conjunction with the CALCOMP Graph Plotter, utili- 


zing the graph-plot program. 
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